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We analyse the thermal evolution of the Universe in the Inert Doublet 
Model for three viable regions of Dark Matter mass: low, medium and high 
DM mass. Those three regions exhibit different behaviour in the possible 
types of evolution. We argue that the quartic self-couplings in IDM are 
I— I significant parameters for the astrophysical analysis. 
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, 1. Thermal evolution of the Universe in IDM 

^ The Inert Doublet Model (IDM) [H [2] is a Z2-symmetric 2IIDM, which 

I> for a special set of parameters may provide the Dark Matter (DM) candi- 

date. The model contains two scalar SU{2) doublets: a "standard" scalar 
(Higgs) doublet and a "dark" scalar doublet <I>5 is responsible 

for the electroweak symmetry breaking and masses of fermions and gauge 
bosons as in the Standard Model (SM), while does not receive vacuum 

^— I expectation value (v.e.v.) and does not couple to fermions. In the model 

the discrete D-symmetry of the Z2 type is present: 

> D: ^s'^^s, ^D^-'^D, fields A 5M fields. (1) 

X 

^ All the degrees of freedom of the dark doublet are realized as the 

massive D-scalars: two charged and two neutral Dh and Da- They 
possess a conserved multiplicative quantum number, the odd D-parity, and 
therefore the lightest particle among them can be considered as a candidate 
for the DM particle. 
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The D-symmetric potential V, which can describe IDM, is: 



2 



2 

$^<I>d) +/i.c. 



, (2) 



with ah parameters real and A5 < [3j. Positivity conditions imposed on 
the potential guarantee that the extremum with the lowest energy will be 
the global minimum of the potential (vacuum). Relevant conditions are: 
Ai,2 >0, R + 1>0; R = A345/\/AlA^ , A345 = A3 + A4 + A5. 

The Yukawa interaction of SM fermions ■0 / with only one scalar doublet 
have the same form as in the SM with the change $ — t- <I>5 (Model I for 
a general 2HDM). 

We consider thermal evolution of the Lagrangian, following the approach 
presented in [H El [3] . In the first approximation the Yukawa couplings and 
the quartic coefficients of V are constant, while the quadratic parameters 
mf- {i = 1, 2) vary with temperature T as follows [31 [6]: 



mi{T) = mi-c,T\ 

„ _ 3A1+2A3+A4 I 3g^+g'2 , _ 3A2+2A3+A4 I 3g^+g'- 

g -rg-r2'''2 6 """8 



(3) 



In this work we limit ourselves to positive ci, C2 as we consider only the 
restoration of EW symmetry for high T (the negative values of m^]^(r) and 
771,22(7") for high enough T). 

As the Universe is cooling down the potential V Q, with temperature 
dependent quadratic coefficients ([s]), may have different ground states |3j. 
The general form of the neutral extremum is: 

( °, ) • ("^-i+^y- w 

EW symmetric extremum, [EWs^ is realized if f£) = tig = 0. Here all 
fermions and bosons are massless and EW symmetry is conserved. 

Inert extremum 1\ can be realized if vn = 0, w| = = m^^^/Ai. If l\ 
is a vacuum then in the scalar sector there are four dark scalar particles 
Djy, Da,D^ and the SM-like Higgs particle hs- The lightest dark particle 
is stable and so it is a good DM candidate. Assuming that DM particles 
are neutral, we consider such variant of IDM in which Dh is a DM candi- 
date, meaning M£)±, Md^ > Md^^. Various theoretical and experimental 
constraints apply for the IDM (see e.g. [3 El |9l EOl H [HI ES]). EWPT 
and collider data (LEP II, Tevatron, LHC) constrain the allowed regions 
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of the masses of hs and dark scalars. The rehc density measurements and 
the direct detection experiments can be used to constrain the DM mass and 
the DM-Higgs self-couphng A345. However, they don't constrain the DM 
quartic self-couphng A2. 

Inertlike extremum I2 is mirror-symmetric to /i as vs = 0, vjj = v'^ = 
771,22 /A2. Fermions are massless at the tree-level (Model I), while gauge 
bosons are massive. There are four scalars Sh,Sa,S^ (no DM candidate 
as D-symmetry is spontaneously violated) and the Higgs particle with 
no interaction with fermions. 

Mixed extremum M is a standard 2HDM extremum with vs,vd 7^ 0. 
Fermions and bosons are massive and there are 5 Higgs particles: CP-even 
h and H, CP-odd A and charged H^, none of them can be a DM candidate. 

We assume that today inert phase Ii with the DM candidate Dh is real- 
ized. However, the sequences of transitions between different vacua (called 
here rays) were possible in the past [3]. 

There are three types of sequences that start in the EWs symmetric 
phase and end in the inert phase. First one is a single phase transition 
EWs — )• /i. It is realized by rays Ia,b,c (for R>1, 0<R<1 and 
—1<R<1, respectively), rays na,b (i? > 1 and < i? < 1) and ray HI 
(only for R > 1). The difference between those rays is the status of the I2 
extremum: for ray I it's not an extremum; for ray H it is an extremum, but 
not a minimum; for ray HI it is a local minimum, but not the global one. 

Second type of sequence, EWs — I2 — )■ Ii, can be realized only if i? > 1 
and is represented by the rays IV and V. In this case the EWSB is a phase 
transition of a 2nd-order, while the last transition I2 — >■ Ii is of the Ist-order. 

Only for R > 1 there is an unique opportunity of coexistence of minima 
(vacuum Ii and local metastable minimum I2) for rays HI, IV and V. For 
ray IV the coexistence is temporary and the local minimum I2 disappears 
for low temperatures, while for rays HI and V it still exists for T = 0. 

In the < -R < 1 case there is a possibility of having three phase 
transitionss in the sequence EWs — )• I2 — )• M — )• /i (ray VI). All transitions 
here are of the 2nd-order. 

2. Phenomenological analysis 

For phenomenological studies it is useful to chose the physical masses 
Mhg, Mdjj, Md^, Md± and the scalar self-couplings as an input parame- 
ters. In this analysis we chose two self- couplings, A345 and A2, which have 
different properties and play different role in the analysis [14]. 

A345 is a triple and quartic coupling of the DM partricle and SM-like 
Higgs, i.e. DnDnhs or DuDnhshs- In wide range of DM mass this param- 
eter governs the main annihilation channel into pair of fermions via Higgs 
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,l345 A345 A345 

(a) low DM mass (b) medium DM mass (c) high DM mass 

Fig. 1: Possible rays for different regions of DM mass. Vertical bounds 
denote region allowed by WMAP measurements; region A is exluded by 
positivity constraints, in region B Ii is only a local minimum. 



exchange: DhDh hs ^ ff with the cross-section a cx \^^^/{4:Mf)^ — 
M|^)^. For this reason this parameter, along with the DM mass, influ- 
ences strongly the value of the DM relic density QoAih'^- It also plays an 
important role in the direct detections experiments, as DM-nucleon elastic 
scattering cross-section is given by (Jdm.n oc X\^^/ {Mj^jj + M^)'^ [2j. 

The remaining self-coupling, A2, is a quartic coupling of DM particle. 
For this reason the exact value of A2 does not influence ri^Af^^ directly. 
However, this parameter limits A345 through the positivity constraints and 
is important for the type of evolution. 

Below we present the analysis done in the (A345, A2) phase space for the 
three regions of DM mass for chosen benchmark points: 

1. = 5 GeV,MB^ = 105 GeV,M^± = 110 GeV, M;,^ = 120 GeV, 

2. Mdh = 50 GeV, Mdj, = 120 GeV, M£,± = 120 GeV, Mhs = 120 GeV, 

3. Mdh = 800 GeV, Md^ = 801 GeV, M^i = 801 GeV, Mhg = 120 GeV. 



Figures la|lb|lc show the possible rays for each region, as well as the 3(T 



WMAP-allowed regions 0.085 < l^DAf/i^ < 0.139 [15]. Note that the region 
A is excluded by the positivity constraints and in the region B Ii is only 
a local minimum and not the vacuum. Each ray is realized in the separate 
region of (A345, A2) phase space. Furthermore, different types of evolution 
are possible in the cases of low, medium and high DM mass. 

2.1. Low DM mass 



The low mass region, which resembles the singlet scalar DM, has been 
shown to fit into the CoGeNT, DAMA/Libra and CRESST-II signal [HKITl 
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[ig m [20], however it appears to be excluded by the XENONIOO results 
|21j . In this region the DM particle is much lighter than all other scalar 
particles with (4- 8) GeV and ~ ^^d± ~ 100 GeV [8]. Large 

mass splittings between the Dh and other scalar particles do not allow for 
the coannihilation. To have the correct WIMP cross-section and proper 
relic density rather large A345 is needed [H]. 

In this region the possible types of evolution are limited to three rays 
only (ray Ic, lib, VI, figure la) and there is no coexistence of minima. 

Notice, that to fit into the WMAP data we need not only large A345, 
but also large A2 ~ 1. The smaller values of A2 are excluded by positiv- 
ity constraints (region A) or I2 vacuum (region B). However, larger A345 
corresponds to the lower temperature of the final phase transition. In this 
example for the sequence EWs I2 ^ M ^ Ii it occurs at Tm^Ii = 6 
GeV, so Mjj^ Tm^Ii ■ The recent analysis [6] shows that in this case the 
lowest order of the thermal corrections to V is not sufficient. 



2.2. Medium DM mass 

In this region the DM mass is of the order of ~ (45 — 160) GeV. 

Mass splitting between Dh and D^ should be large: M£,± —M^^ w (50—90) 
GeV. Constraints for Mq^ —Mq^ have been derieved in p.lj . If this value is 
large with Mq^ ~ -^d± then there is no coannihilation. For Mo^ — Md^ < 
8 GeV this effect influences strongly the value of DM relic density [51 [9] . 

For medium DM mass, regardless of the exact values of the mass split- 



tings, all rays are possible (figure lb). This is the only case when one can 



have the Ist-order phase transition for rays IV and V - those rays are not 
possible for low or high DM mass. 

In this region Q.^Mh'^ is very sensitive to the exact value of M^^ and 
mass splittings. Therefore, we cannot make a general statement that a 
certain ray will always give a proper relic density [H]. However, some prop- 
erties of this mass region are independent on M^)^. For example, complex 
sequences (rays IV- VI) require rather large A345. This however leads to the 
similar problem as in the low DM mass case: the temperature of the final 
phase transition is lower and especially for ray V further thermal corrections 
to the potential are needed. 



2.3. High DM mass 

As shown in [9] in the high mass region the mass of the dark matter 
particle should be over 500 GeV, with the small mass splittings between the 
dark particles. Also the perturbative unitarity may give the relevant con- 
straints [22]. To have a proper relict density we need rather large absolute 
values of IA345I and A2 [T^. Only three rays are possible; they correspond 
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Fig. 2: Temperature evolution of masses of scalar particles in different se- 
quences. Mdh = 50 GeV, = 120 GeV, M^± = 120 GeV, Mhg = 
120 GeV, A345 = 0.1945. Notation: D = D^, Da, S = S^, Sa- 



to the sequence with a single phase transition EWs — )■ /i and differ only by 
Ri 
20) 



the value of R (rays la, lb, Ic, figure Ic). Other types of evolution require 
A ^ 0(10 



2.4- Role of X2 self-coupling 

X2 self-coupling is a significant parameter in IDM - as shown in the 
previous sections. It not only limits A345 through the positivity constraints, 
but also is important for the type of evolution. In this section we fix the 
scalar masses and A345 self-coupling and let A2 vary. Depending on the value 
of A2, different types of evolution are realized. 

Figures 2apb show the thermal evolution of the mass parameters dur- 



ing evolution of the Universe [23j. In case of ray V (figure 2a) there are 
two phase transitions. EWSB into I2 phase happens for T = 134.8 GeV. 
Dashed line shows the apperance of the local minimum Ii during the inert- 
like phase of the evolution. Ii becomes a global minimum after Ist-order 
phase transtion which takes place for T = 83.7 GeV. From this point I2 is 
a local minimum that still exists for T = 0. 



Figure 2b shows the evolution according to ray VI. Here, after EWSB at 
T = 126.7 GeV Universe enters the inertlike phase with massless fermions 
and massive gauge bosons. This minimum becomes a saddle point for T = 
119.4 GeV and the 2nd-order transition to the M phase takes place. This 
phase is short-lived and soon, at T = 119.0 GeV, there is another 2nd-order 
transition into Ii phase. There is no coexistence of minima at any point in 
time during evolution. 

Notice, that although the temperature of EWSB is similar in both cases, 
the final phase transition happens at different temperatures and it's much 
lower for ray V. As discussed, this ray is the most likely to require further 
thermal corrections to the potential [6]. 
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3. Conclusions 

In this paper we studied the temperature evolution of the Universe in 
IDM. We also discuss the significance of the quartic self-coupling A2- This 
parameter does not influence DM relic density directly and it cannot be 
accessible in the colliders. However, it is related to the value of A345 coupling 
through the positivity constraints. It also plays an important role in the 
evolution, as its different values lead to the different types of the evolution 
of the Universe. 
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